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ABSTRACT 

1-loop diagrammatic calculations of cross sections and decay widths of neutral 
Higgs bosons in the minimal supersymmetric standard model are reviewed and 
compared with compact expressions in the effective potential approximation. 

1. Introduction 

In order to experimentally detect possible signals from the Higgs sector of the 
minimal supersymmetric standard model (MSSM), detailed studies for the decay 
and production processes of Higgs bosons are required. As discovered several years 
ago a'iH, radiative corrections in the MSSM Higgs sector are large and have to be 
taken into account for phenomenological studies. Three main approaches have been 
developed to calculate the 1-loop radiative corrections to the MSSM Higgs boson 
masses, production and decay rates: 

(i) The Effective Potential Approach (EPA) I. 

(ii) The method of Renormalization Group Equations (RGE) i. 

(iii) The diagrammatic calculation in the on-shell renormalization scheme (Feynman 
Diagram Calculation, FDC) lii'lii: The masses are calculated from the pole 
positions of the Higgs propagators, and the cross sections are obtained from the 
full set of 1-loop diagrams contributing to the amplitudes. 

The method (iii) is technically involved, but it is the most accurate one at the 1-loop 
level and can be used as a reference frame for simpler approximations. The searches 
for Higgs bosons at LEP and studies for the future searches at higher energies 
conventionally make use of the very compact formulation in the effective potential 
approximation. 

This talk gives an overview on the neutral MSSM Higgs sector at the 1-loop 
level. The results for the cross sections of neutral Higgs production processes in 
e + e _ collisions and for the neutral Higgs decay widths are discussed in a complete 
diagrammatic calculation and compared, where possible, with the corresponding ones 
of the compact EPA approximation. 
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2. One-loop calculations 



The tree level potential for the neutral MSSM Higgs bosons can be written as 
follows: 

V = m\Hl + m\El + t l3 {m\ 2 H\H{ + H.c.) + 9 —^{Hl - tf 2 2 ) 2 + 9 ^{H x H 2 f . (1) 

Diagonalization of the mass matrices for the CP-even and the CP-odd scalars, fol- 
lowing from the potential ([!]), leads to three physical particles: two CP-even Higgs 
bosons H°,h° and one CP-odd Higgs boson A , and defines their tree-level masses 
m hy m Hy m A and the mixing angles a, (3. For a sytematic 1-loop caclulation, the free 
parameters of the Higgs potential ml, m|, m\ 2) g, g' and the two vacua vi, v 2 are 
replaced by renormalized parameters plus counter terms. This transforms the poten- 
tial V into V + SV, where V, expressed in the renormalized parameters, is formally 
identical to ([[]), and 5V is the counter term potential. The counter terms are fixed 
by seven renormalization conditions. In the on-shell scheme they can be chosen as 
follows: 

• the on-shell conditions for My/,z and the electric charge e as in the minimal 
standard model. 

• the on-shell condition for the A boson with the pole mass Ma- 

• the tadpole conditions for vanishing renormalized tadpoles: 

T H + 5t B = 0, T h + 5t h = 

where T B ,h are the sum of the 1-loop tadpole diagrams for H° and h°, and Stn,h 
are the tadpole counter terms following from (1). These conditions ensure that 
Vi, t>2 are the minima of the potential at the 1-loop level. 

• the renormalization of tan/3 in such a way that the relation tan/3 = i^/fi is 
valid for the 1-loop Higgs minima. 

By this set of conditions, the input for the MSSM Higgs sector is fixed by Ma and 
tan j3, together with the standard gauge sector input M w ^ z and e. The last condition 
on tan f3 can only be imposed in connection with an appropriate field renormalization 
of the two Higgs doublet fields. Together with the gauge field renormalization one 
has four extra renormalization constants which can be fixed as in the standard model 
gauge sector 0, extended by two more conditions for the Higgs sector. The latter two 
have been treated in two slightly different ways i'i in the literature; physical results, 
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however, differ only marginally by unobservably small terms. The corresponding 1- 
loop physical Higgs boson masses M^Mfj are obtained as the pole positions of the 
dressed scalar propagators. 

In the EPA, the tree level potential is improved by adding the 1-loop terms . The 
1-loop potential V^ 1 ) is rediagonalized yielding the 1-loop corrected physical masses 
Mn,Mh and the effective mixing angle a e // • These improved masses and a e // are 
used in the Born formulae for production and decay rates of Higgs bosons. In the 
approximation keeping only the dominating terms ~ mf, the expressions read i: 
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The approximate effective mixing angle a e /j is determined by 

-{M\ + M 2 Z ) sin /3 cos/? + A t 
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These formulae contain the masses mi of the top squarks, the Higgs mixing param- 
eter /i of the superpotential, and the non-diagonal entry A t in the stop mass matrix. 
This matrix is diagonal for A t + fi cot (3 = 0. In this special case we have at = Xt = 0, 
and only the e t term contributes. 

Recently the leading 2-loop corrections to the CP-even MSSM Higgs boson masses 
have been investigated, based on the EPA and RGE methods 0. The main conclusion 
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is that 2-loop corrections are also significant and tend to compensate partially the 
effects of the f-loop corrections. 

3. Production cross sections for e + e~ — > Zh°(H°), A°h°(H°) 

In this section the results for Z H°(Z°h ) and A H°(A°h ) production are shown 
as derived from the complete 1-loop FDC, and the quality of the corresponding EPA 
results is discussed. The formulae for the cross sections obtained in the FDC differ 
from the Born expressions not only by the corrections to the masses and to the angle 
a, but also by new form factors and momentum dependent effects (see @@ for analytic 
expressions) . 

For the calculations of the cross sections we need the full set of 2-, 3- and 4-point 
functions. In Fig. |I] the diagrams contributing to e + e~ — > Z°h°(H°) are collected. 
The diagrams for e + e~ — > A°h°(H°) can be obtained by changing Z° into A on the 
external line and skipping the diagrams i), j). 




Figure 1: Classes of diagrams contributing to the e + e — ► Z°h°(H°) process in the FDC 
approach. 

In the figures of this section the set of parameters (in GeV): Ma = 200, M2 = 
1000, /i = 500, M sl = 300, M sq = 1000, A t = A b = 1000 is used as an example. // is 
the parameter describing the Higgs doublet mixing in the MSSM superpotential. M2 
denotes the SU(2) gaugino mass parameter. For the U(l) gaugino mass we use the 
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value Mi = | tan 2 9wM 2 , suggested by GUT constraints. M sq , M s i, A t and Ab are the 
parameters entering the sfermion mass matrices. For simplicity we assume a common 
value M sq for all generations of squarks, and a common M s i for sleptons. 




Figure 2: Comparison of the cross sections a(e + e~ -> Z°h ,A°h ) obtained in the EPA 
and FDC. ^Js = 500 GeV. 

In the conventional parametrization, Fig. ^ shows the production cross 

sections <r(e + e~ — ► Z°h°, A°h°) at y/s = 500 GeV. For the choosen set of parameters 
the numerical differences can reach 30% at y/s = 500 GeV. Note, however, that in 
the region of large cross sections the EPA accuracy is better (20% at 500 GeV). The 
situation of H° production is very similar. The differences between EPA and FDC 
become more important with increasing energies, exceeding 40% at 1 TeV. Also the 
effect of the additional form factors in the FDC grows, which modify the angular 
dependence of the cross section compared to the effective Born approximation. More 
detailed discussions can be found in ref. §0. 

A more physical parametrization of the cross sections is given in terms of the 
two Higgs boson masses Ma and (or M# ), instead of the formal quantity tan/?. 
This parametrization is more involved in the calculations, but it has the advantage 
of physically well defined input quantities avoiding possible confusions from different 
renormalization schemes. Varying e.g. Mh {Ma and other input quantities fixed) we 
obtain tan/3 and <Jzh, &AH as functions of M H . Significant differences can occur for 
the cross sections, as displayed in Fig. [3] where the predictions of EPA and FDC for 
the gzh and o ah are plotted as functions of Mh ■ The typical size of the differences 
between the methods is 10-20% for ^Js = 500 GeV, but they may became quite large 
(60%) for the process tx(e + e~ — > Z°H°). In other cases, they are of the order 10-20%. 
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Figure 3: Comparison of the cross sections <r(e + e~ — ► Z H°,A°H ) versus Mh in the 
EPA and FDC. = 500 GeV. 

The variation of the SUSY parameters: sfermion and gaugino masses, /x parameter 
and sfermion mixing parameters, does not have a large effect on the size of the 
differences between the EPA and FDC. Hence, the figures represent typical examples. 

Summarizing this section, comparisons between the FDC and EPA predictions 
have shown that at ^/s = 500 GeV the EPA has an accuracy of typically 10-20% in the 
parameter regions where the cross sections are large. The differences become larger 
with increasing energy, where also modifications of the Born-like angular distributions 
are more visible. The use of the physical input variables Ma, Mh or Ma, Mh avoids 
ambiguities from the definition of tan j3 in higher order, but the observed differences 
remain of the same size. For a better accuracy, the full FDC would be required. 

So far the leading 2-loop terms have not been incorporated. They would im- 
prove the 1-loop FDC results in the same way as the approximations and thus do 
not influence the remaining differences which can only be obtained by an explicit 
diagrammatic calculation. 

Loop-induced pair production e + e~ — > h h ,H H ,h H ,A°A of neutral Higgs 
bosons have also been studied recently 0, as well as associated Higgs-photon produc- 
tion e + e~ — > 7/i°(if °, A ) S The general result is that the cross sections in the MSSM 
are not enhanced by the extra non-standard particles in the loops, a(MSSM) < 
a(SM). In the decoupling limit, for heavy SUSY particles, the standard model results 
are recovered. For detailed studies of the virtual SUSY effects in the loop-induced 
Higgs-7z7 and Higgs-77 couplings see ref. 0. 
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4. Decays of neutral Higgs bosons 



The decay widths (the branching ratios, respectively) as well as the mass-width 
correlations are quantities which can help to differentiate between Higgs bosons of 
different origin. Except for a small region of the parameter space, the light neutral 
Higgs of the MSSM decays predominantly into 6-quarks and r-leptons; the heavier 
ones H° and A can have significant decay modes also into top quarks, scalar quarks, 
and neutralinos/charginos. In a certain region of the parameter space, also the decay 
h° — > A A is allowed. Loqrj-mediated decay processes are the hadronic decay modes 
into gluons and gluinos 
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4-1. Fermionic decays 

For the important fermionic decays both electroweak and QCD corrections 000 
have been calculated. The standard QCD corrections for <p — > qq, </> = h°, H°, A , 
are large. The bulk can be absorbed into the running quark mass by replacing the 
pole mass according to m q — ► m q (M^). The SUSY-QCD corrections arising from 
virtual gluinos and squarks can also become remarkably large, in particular for 
large values of fi and tan/? where they can reach up to 30%. 

The set of 1-loop electroweak corrections to h° — > // can be summarized in terms 
the following decay amplitude: 

A ( k m = ^ (r» - Mi _ ^§Lm r «) < 5 > 

with the renormalized self-energies £ and 3-point vertex functions Th,u- The am- 
plitude for H° — > // is obtained by interchanging h «-> H. The wave function 
renormalization Z^h) is the finite residue of the h° (H°) propagator. The amplitude 
for A — ► // is given by the renormalized vertex Ta alone, due to the renormalization 
condition Za = 1 i. For the absolute decay widths, in order to have the correct 
normalization in terms of the Fermi constant Gf, the MSSM correction to the muon 
lifetime, i.e. the quantity Ar, has to be taken into account 0. It drops out in the 
branching ratios. The inclusion of the mixing term in eq. ([5]) corresponds essentially 
to the rediagonalization of the mass matrix in the EPA. In the EPA, one obtains the 
improved decay amplitude by using the EPA masses and the effective mixing angle 
a e ff , eq. (f|), in the Born expression for the vertex T h , with Z h = 1 and T, hH = 
in eq. (||). As shown in I, the EPA is a very good approximation of the full 1-loop 
result for the fermionic branching ratios, with exception of extremely low values for 
tan/5 (see Fig. [§. 
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h° -> ff branching ratios 
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Figure 4: Fermionic branching ratios Higgs decays into fermion pairs (from 0). EPA: dashed 
lines; complete 1-loop: full lines. M denotes the SU(2) gaugino mass. 
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4-2. The decay h° -> A°A° 



A specific consequence of the large radiative corrections to the h° mass is the 
possibility of having > 2Ma, which makes the decay h° — > A°A° kinematically 
allowed at 1-loop order for low values of tan j3 and Ma- In the allowed region it turns 
out to be the dominant decay mode, with branching ratios of 0.8-0.9. The decay 
width was obtained in the EPA (last reference of @), by use of the RGE il, and by a 
complete diagrammatic 1-loop calculation il. A simple approximate formula for the 
decay amplitude, which reproduces the full 1-loop decay width at an accuracy within 
typically 10%, is given by Ej 
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which arises from the vertex correction diagrams with virtual top and stop quarks. 
It contributes significantly to the decay amplitude in particular for tan/3 ~ 1, where 
the improved Born approximation [the first term in eq. (^|)] is very small. The second 
term in AV^ is also sizeable for intermediate masses of the top squarks. For details 
see ref. 0. 



4-3. Higgs decays into s quark pairs 

In a large part of the MSSM parameter space, the decays into squark pairs 
H°, A — > qq can be the dominant decay modes @. QCD corrections from gluons 
and gluinos have been calculated recently HI. They are sizeable (up to 50%) and 
should be taken into account for phenomenological studies. 
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